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Abstract 

Based on signaling process on complex networks, a method for identification com- 
munity structure is proposed. For a network with n nodes, every node is assumed to be 
a system which can send, receive, and record signals. Each node is taken as the initial 
signal source once to inspire the whole network by exciting its neighbors and then the 
source node is endowed a nd vector which recording the effects of signaling process. So 
by this process, the topological relationship of nodes on networks could be transferred 
into the geometrical structure of vectors in nd Euclidian space. Then the best partition 
of groups is determined by F-statistic and the final community structure is given by 
Fuzzy C-means clustering method (FCM). This method can detect community struc- 
ture both in unweighted and weighted networks without any extra parameters. It has 
been applied to ad hoc networks and some real networks including Zachary Karate 
Club network and football team network. The results are compared with that of other 
approaches and the evidence indicates that the algorithm based on signaling process is 
effective. 

Keyword: Complex network. Community Structure, Signaling Algorithm, FCM 
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1 Introduction 



The study of complex networks has received an enormous amount of attention [TJ from 
the scientific community in recent years. Physicists in particular have become interested 
in the study of networks describing the topologies of wide variety of systems, such as the 
world wide web[4j, social and communication networks O [6], biochemical networks and many 
more[7]- One such problem is the analysis of community structure found in many networks. 
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Distinct communities or modules within networks can loosely be defined as subsets of nodes 
which are more densely linked, when compared to the rest of the network^ [9]. Such 
communities have been observed, using some of the methods, in many different contexts, 
including metabolic networks, banking networks and most notably social networks. As 
a result, the problem of identification of communities has been the focus of many recent 
efforts. Community detection in large networks is potentially very useful. Nodes belonging 
to a tight-knit community are more than likely to have other properties in common. Besides, 
these communities may probably be functional groups, which provide us valuable reference 
to our study in many other fields. In recent studies, the scientists have designed many 
different algorithms [i [i[ini[IIl[Il[Il[Il[I3[ia[ni[ISl[ini[ini[lIl[ll[23[21] (see [12] 
as a review) to detect the community structures. These algorithms can be divided into 
categories. Some algorithms are designed according to maximal modularity Q. Some are 
designed based on topology structures (betweenness, degree, or clustering coefficient). The 
last is designed according to the dynamical properties of the network. 

Communities within networks can be defined as subsets of nodes which are more densely 
linked, when compared to the rest of the network. Modularity Q is an index advanced 
by Newman and Grivan[25' as a measurement for the community structure. It gives a 
clear and precise definition of characteristics of the acknowledged community and have very 
successful application in practice. So it leads to many other algorithms brought forward 
to divide a community by maximizing Modularity Q. Unfortunately, maximizing Mod- 
ularity Q has been proven a NPC problem [26j. which makes it unable to work out the 
partition corresponding to maximal Q in a network which has lots nodes. Actually many al- 
gorithms for maximizing Q are usually heuristic algorithm. Besides, Modularity Q has been 
proven strictly that, as an index to measure the community structure, it tends to combine 
the little communities rather than identify them successfully in the networks with definite 
communities [57]. Though Modularity Q is proven to have the above-mentioned inherent 
defects, it is still the successful index to measure a network for the moment. Therefore, lots 
of work for detecting communities are dependent on the index Q. 

Among the algorithms based on network topology, we want to introduce briefly the spec- 
tral analysis method and GN algorithmjS]. GN algorithm was proposed by Grivan and 
Newman. It gives the division of network by remove links. The links with largest between- 
ness are removed one by one in order to split hierarchically the graph into community. But 
GN algorithm alone can only give the dendrogram concerning the network structure finally. 
It could not give the best partition directly. If we want to get the best partition we have to 
depend on Modularity Q or other indexes to work it out. While the principle of the spectral 
analysis method [SIT is based on the theory of eigenvector of matrix. When a network is 
partitioned into two communities with the pre-fixed sizes, the best partition is to make the 
number of edges between the two communities being minimum. In fact, this problem is 
closely related to the eigenvector corresponding with the Fiedler eigenvector (the eigenvec- 
tor correspond to second minimum eigenvalue) of Laplacian matrix of a network. Relatively 
speaking, spectral analysis method is the most math theory-based approach. But it's still 
unable to get the best partition. Besides, it requires to know the sizes of the two sub- network 
beforehand. Although some methods are proposed to solve the two defects, like getting the 
best partition by using Q, ascertaining the sizes of the two sub-network by the sign of the 
elements of Fielder eigenvector and so on, the inherent defects of spectral analysis method 
have not been solved perfectly and completely yet. 
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There are still other algorithms based on the dynamics on networks, among which random 
walks method and circuit approach methodrTT] will be briefly discussed here. In the 
random walks method, each node contains a walker initially. Then each walker will randomly 
choose a neighbor of the node it currently stand on to localize. This is a Markov process. 
After a period of time, the possibility will be higher that the walker reach another node 
belonging to the same community of the node the walker stood on at the beginning. And 
this possibility can be directly regarded as the possibility of pair of nodes in the same 
community and through it, a dendrogram can be got, then partition can be made by the 
aid of Modularity Q. 

When using random walk method to detect communities, it's difficult to specify the 
optimum random-walking time. And the best partition is dependent on the external index 
Q usually. The principle of circuit approach method is to regard the edges of the network 
as the resistances, and add voltage to the adequate nodes of the network, then work out the 
voltage of each node by Kirchhoff's law. Nodes with the similar voltages are regarded to 
exist in the same community more probably. At the same time, define the external indexes 
such as tolerance to realize the partition of the network. 

In this paper, we propose a new algorithm for identification communities based on sig- 
naling process on network. In this approach, every node is viewed as a system which can 
be inspired. It can send, receive, and record signals. In the initial, a node is selected as 
the source of signal. We give it an initial unit signal and other nodes with the signal of 
zero. Then the source node send the signal to its neighbors and itself first. Afterwards, 
the nodes with signals can also send signals to their neighbors and themselves. One thing 
should be mentioned in this signaling process is that the node can record the amount of 
signals it received, and at every time step, each node sends its present-owning signals to its 
neighbors and itself. After the inspiration of a certain T time steps, the signal distribution 
over the nodes could be viewed as the influence of the source node to the whole network. 
For a network with n nodes, signal distribution can be characterized by a nd vector. 

If a network has n nodes, we can get the influence of every node by the same operation. 
The results are given by n nd vectors. Generally speaking, the source node should influence 
its community first then through the community to influence the whole network. So natu- 
rally, compared with the nodes in other communities, nodes of the same community have 
similar influence toward the whole network. And the difference of influence could be given 
by the nd vectors. 

Thus, by the above signaling process on networks, the topological structure of nodes 
is converted into the geometrical relationships of vectors in nd Euclidian space. We can 
get the community structures of nodes by clustering these nd vectors. Actually, there are 
already a lots of methods to cluster vectors in Euclidian space. Here we chose Fuzzy C- 
means clustering method (FCM) assisted by F statistic[ST] to get the best partition of the 
communities. F statistic is developed in mathematical statistics. It describes the best 
partition as the one with the shorter average distance between the vectors inside the same 
community, and the larger distance between vectors of different communities. After getting 
the best number of groups by F statistic, we can work out the communities by FCM. With 
the aid of F statistic, the method presented in this paper can detect community structures 
in complex network without any extra parameters. 

Some problems related to above method are also discussed, including the optimum time 
steps T of inspiration and the generalization of the method to weighted networks. Then 
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we applied the method to detect the communities in ad hoc and some real networks. Its 
precision and accuracy are obtained and compared with some other algorithms. The results 
indicate that the method based on signaling process performs good. 

2 Method Based on Signaling Process 
2.1 Basic Algorithm 

A. Signaling Process. For a network with n nodes, every node is assumed to be a system 
which can send, receive, and record signals. One node can only affect its neighbors which 
will affect their neighbors too in the same way. Finally, each node will affect the whole 
network. In general, one node will affect its community first and then the whole network 
via its community. So we can safely conclude that the nodes in the same community will 
affect the whole network in a similar way. 

At the beginning, we select a node as source and let it has one unit of signal and the 
other nodes have no signal. Then let the source node send signal to all of its neighbors and 
itself. After the first step the node and all its neighbors have a signal. In the second step, all 
the nodes which have signal will send it to their neighbors and themselves. Every node will 
record the amount of signals it received and then it will send the same quantity of signals 
in the next time step. In this way, the process will be repeated constantly on the network. 
After T time steps, we can get a nd vector that records each node's signal quantity which 
represents the effect of the source node. The signaling process is sketched out in FiglT] by 
simple network with 5 nodes. Choosing every node as the source node respectively, we can 
get n such vectors. The purpose that we let each node sends signal or signals to itself is to 
take account of the historical effects. This has been proved to be helpful to distinguish the 
amounts of signals between the nodes in the community and outside in a relatively short time 
period. Standardizing the n vectors, then the distance of each pair of vectors will represent 
the similarity of the corresponding nodes. Using this kind of similarity the communities can 
be detected. 

Actually, the above signaling process could be described by a simple but clear math- 
ematical mechanism. Suppose we have a network with n nodes, it can be represented 
mathematically by an adjacency matrix A with elements Aij if there is an edge from i to 
i and otherwise. Then the column i of matrix V = (I + A)-^ will represent the effect of 
source node i to the whole network in T steps. In order to get the relative effect, we should 
standardize each column of matrix V. Assume the column i of V is I^ = (f a, i'i2---'f^m), 
then the Vi can be standardized as Ui = {un, Ui2...Uin), here Uij — ^S^. Then to partition 
the network which has n nodes is equivalent to cluster n vectors Ui,U2, ■ ■ ■ ,Un in Euclidian 
space. 

B. Fuzzy C-mean clustering. It is well known that there are many clustering methods 
and algorithms for the vectors in Euclidian space. In this paper, we choose the inexpensive 
fuzzy C-mean clustering algorithm (FCM)|3T] to detect communities for the vectors given 
by signaling process. FCM is described as following. 

1. Set C as the number of communities to partition. 
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Figure 1: Sketch of signaling process, (a) We set a signal on node 1 and other nodes have 
no signal initially. (b)In the first step node 1 sends one signal to all its neighbors which 
are nodes 2,3,4 and 5 and itself. Then all nodes have 1 signal. (c)Next, each of them will 
send one signal to their neighbors and themselves respectively at the same time. After the 
second step, node 1 has 5 signals, both node 2 and 3 have 4 signals and node 4 and 5 have 
3 signals. Vector [5,4,4,3,3] represents the effect of the node 1 to the whole network in 2 
steps. (d)Then every nodes send the same amount of signals as they received in the last 
step to their neighbors and themselves. 

2. Randomly choose C vectors for the C communities as their barycenters. 

3. Randomly choose a vector. The vector will belong to the community when the distance 
between the vector and the barycenter of the community is minimum among all the 
barycenter of communities. 

4. Re-compute the communities' barycenters which have added a vector or deleted a 
vector. 

5. Repeat step 3 to step 4 until all the barycenters cannot be modified. 

We know that there are many definitions of distance. In our algorithm we choose the 
normal definition-Euclidian distance to measure the similarity between vectors of nodes. 

C. F Statistic. At the first step of fuzzy C-mean clustering algorithm, we must set 
an extra parameter C which presents how many clusters we will partition. Here we use F 
statistic[2I] to estimate the proper C. Now let we have a glance at F statistic. Suppose 
U = {ui, U2, ■ ■ ■ , Un} is the set of vectors of all nodes and Uj = (xji, Xj2, ■ ■ ■ , Xjn), here Xjk 
is the kth character quantity of Uj. Suppose C is the number of communities and rii is the 
number of nodes of ith community. We name all the nodes' vectors of the ith community as 
u\,U2,---, mJj. . Let x], — 5I]J=i k — 1, 2 • • • , ti be the mean characters of ith commu- 
nity, = {x\,xl, • • • , xjj be the ith. community's barycenter and u = {xi )be all 
the nodes' barycenter, here Xk = ^jk = Ij 2, • • ■ , n). Then F statistic is defined 
as 
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Figure 2: F statistic as the function of number of clusters C. Plot shows the changes of F 
statistic when T = 3 and (kinter) changes from 4 to 8. It shows that when (kinter) is smaller 
than 8, F statistic can identify the right number of communities. When the community 
structure is clearer, the maximal value of F statistic is very distinct. 

where || —u \\= \/'^'k=ii^k ~ is the distance between M*and m, and || — -«* || is the 
distance between m^- node of zth and the barycenter . The numerator of F signifies the 
distance of inter-communities and the denominator the distance of intra-communities. So 
the F is bigger when the difference distance of inter-communities is bigger and the difference 
intra-communities is smaller. We can image that when F achieve the maximum we can get 
the best partition. 

We use binary ad hoc networks which contains 128 nodes and 4 groups and proceed 
the signalling process as above to test the F statistic. The results show that F statistic is 
very efficacious. On the weighted ad hoc networks, the results are similar with binary ones. 
When the community structure is clearer the maximal value of F statistic is more distinct. 
The detailed results are shown in Figl2l 

2.2 Some Related Problems 

A. The Most Optimal T. Parameter T is an important factor for the results of community 
identification. We can image that the length T must be sufficiently long to gather enough 
information about the topology of the network and it should not be too long to faint the 
information we have gathered. In order to let majority of nodes can affect the whole network 
and do not to faint the information about the topology of the network, we guess that it may 
be optimal when T is near to the average shortest path of network. In order to demonstrate 
our guess, some numerical experiments are done on binary networks which contains 128 
nodes and 4 groups as above. The results are shown FigO The accuracy of the algorithm 
reach optimum when T is 3 or 4. Of course, we only do some numerical experiments, it is 
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Figure 3: Let the number of communities be 4 always. The plot shows the changes of 
accuracy when T changes form 1 to 10 on binary ad hoc networks. From the plots we can 
see T « 3 is the proper length. And here 3 is near the average shortest paths of each artificial 
networks. 

hard to say the result satisfies all the networks. The random walks method[T4j has also the 
same problem. How to find the most optimal T? We think it is still an open problem. 

B. Time Complexity analysis. The time complexity of computing F statistic is 
O(n^). For a definite C which is the number of communities, the time complexity for Fuzzy 
C-mean clustering is 0{Cn'^). Time complexity of the process of signal diffusion is 0{Tn?) 
when we use the multiplication of matrix to simulate the process. But if we simulate the 
process in a network directly, the corresponding time complexity is 0{T{k + Vjn?), where, 
k is the average degree of node in a network. 

C. Generalization to vifeighted networks. It is easy to generalize our algorithm to 
weighted network. Suppose we have a weighted network with n nodes, it can be represented 
mathematically by an adjacency matrix W with elements Wij . Wij denotes the connection 
strength of node i and j (In some weighted network Wij dose not denote the strength of 
connection, we should transform the weight before the algorithm), then V = (1 + W)-'". The 
rest of the algorithm is same with the algorithm on binary ad hoc network. 

D. Relations with other methods. There two main differences between our method 
and random walks method and circuit approach method. First, we use the signal diffusion 
process to transfer the topology to geometrical structure. The mathematical form is (I+A)-^. 
The distance of each pair of column vectors of the matrix is the intimacy of the corresponding 
pair of nodes. The random walk method gets the intimacy of each pair of nodes by random 
walks. The mathematical form is (diag{-^^ ■ ■ ■ -^)AY' where, di denote the degree of 
node I, diag means the diagonal matrix. Take account of the effect of node degree, it also 
use the Euclidian distance to define the intimacy. The circuit approach method gets the 
intimacy of each pair of nodes by Kirchhoff's law. Adding pressure on the proper two 
nodes, by Kirchhoff's law, we can get the pressure of each node. More close of the pressure 
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of two nodes are, more intimate the two nodes are. Suppose add pressure on nodes 1 
and 2, then pi = 1,P2 — 0, where pi denotes the pressure of node i. The mathematical 
form of Kirchhoff's law is B = {Adiag{-j^^ ■ ■ ■ ^)), (P3,P4, • • ■ ,Pn) — {I ^ B)^^C, where 
B denotes the matrix B with deleting the first and the second columns and rows, C = 
(^^' ^T' ' ' ' ' ^^-^ ■ ^6<^ond, as to the method of clustering, we use the F statistic and 
classical FCM method to partition the vectors. When the F statistic achieve its maximum, 
we get the best partition. The random walks method and the circuit approach method 
are all need the help of other indexes to get the best partitions. One is modularity Q, the 
other is tolerance. So we could say that F statistic and FCM method are all based on the 
geometrical structure of the vector space, but the other two methods need the help of extra 
parameters. Because the random walks method and the circuit approach method both need 
some extra parameters which are very important to the results, so in this paper we will not 
compare our results with that of these two algorithms in different networks. Instead, we will 
compare the accuracy and precision with other famous algorithms which do not require any 
extra parameters. 

3 Results and Comparison with other Algorithms 

In order to investigate our algorithm, the accuracy and precision of our algorithm will 
be compared with Potts algorithm(Potts)T^, Girvan-Newman algorithm (GN)[5S] and ex- 
tremal optimal algorithm(EO)T3j. All these algorithms can be generalized to weighted 
networks |30j. Here we abbreviate ON weighted version as WON and EO as WEO. The 
accuracy and precision are defined in f30' . Accuracy means the consistence when the com- 
munity structure from algorithm is compared with the presumed communities, and precision 
is the consistence among the community structures from different runs of an algorithm on 
the same network. The algorithm's accuracy and precision are calculated by S function [30j 
in this paper. In the following numerical investigations on ad hoc networks, we first get 20 
realizations of artificial community networks under the same conditions. Then we run each 
algorithm to find communities in each network 10 times. Based on these results, using the 
similarity function S*, comparing each pair of these 10 community structures and averaging 
over the 20 networks (average of totally 6*2° x 20 = 900 results) we could get the precision 
of the algorithm. Comparing each divided groups with the presumed structures, we can get 
the accuracy of the algorithm by averaging these 10 x 20 = 200 results. 

3.1 Results on ad hoc networks 

A. Binary ad hoc networks. In order to compare our algorithm with others we first 
test it on computer-generated random graphs with a well-known predetermined community 
structure pS]. Each graph has iV = 128 nodes divided into 4 communities of 32 nodes 
each. Edges between two nodes are introduced with different probabilities depending on 
whether the two nodes belong to the same group or not: every node has (kintra) links on 
average to its fellows in the same community, and {kinter) links to the outer- world, keeping 
{kintra) + (kinter) = 16. Fig[5] shows the rcsults. The precision of our algorithm is better 
than EO and almost the same as Potts and GN. While the accuracy of our algorithm is 
better than GN and almost the same as EO and Potts. 
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Figure 4: Algorithm performance as applied to ad hoc networks with n = 128 and four 
communities of 32 nodes each. Total average degree is fixed to 16. We can see that the 
precision of our algorithm is better than EO and Potts and equal GN and the accuracy is 
better than GN and almost the same as EO and Potts. 

B. Weighted ad hoc networks. In weighted networks, we use similarity link weight 
to describe the closeness of relations between nodes. Under the basic construction of ad 
hoc network described above, the intragroup link weight is assigned as Wintra, while the 
intergroup link weight is assigned as Winter- Similarly with {kintra) + {Winter) = 16, we 
require the link weight on intra and inter links follow the constraint: Wintra + Winter = 2, 
where Wintra ( Winter) IS the average of all intragroup (intergroup) link weights. Here for 
simplicity, we assign the same weight Winter = w to all intergroup links, and assign the same 
weight Wintra = 2 — w to all intragroup links. From Fig[5l we can find that the precision 
of our algorithm is better than WEO and Potts and equal to WGN, the accuracy of our 
algorithm is better than WGN but almost equals WEO and Potts. Even for the case with 
{kinter) < (kintra) but (winter) > (wintra) , Or with Uniform distribution of link weights, we 
can get similar conclusions. 

C. Complete weighted networks. An extreme idealized example is the complete net- 
work. In complete networks, we use uniform distribution of link weights. Weights are taken 
randomly from the interval [{w intra) -0.25, (wintra) +0-25] and [{winter) -0.25, (winter) +0-25] 
respectively, for intragroup connections and intergroup connections. The precision of our 
algorithm is better than WEO and Potts and equal to WGN when (winter) ^ (wintra), 
but its accuracy almost declines to zero when (winter) is greater than 0.9. Fig|S] shows the 
results in detail. 



3.2 Empirical Results on Some Real Networks 

A. Zachary's karate club. Zachary karate club network[28j has bee considered as a simple 
sample for community detecting methodologies [101 [El ED ISl |23l |25]. This network was 
constructed with the data collected observing 34 members of a karate club over a period 
of 2 years and considering friendship between members. We let T = 3 and get the best 
partition which perfectly corresponds to the actual division of the club FiglT] 
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Figure 5: The influence of weight on community structure when (kintra) — (kinter) — 8, 
where Winter changes from 0.05 to 1. We can see the precision of our algorithm is better 
than WEO and Potts and equal to WGN and the accuracy of our algorithm is better than 
WGN but almost equals WEO and Potts. 
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Figure 6: Precision and accuracy of our algorithm and WGN, Potts, WEO algorithms in 
complete networks with presumed communities. When {winter) ^ {wintra), the precision 
of our algorithm is better than WEO and Potts and equal to GN and the accuracy of our 
algorithm is better than GN always but sharply declines to near 0.2 when {winter) is greater 
than 0.9. 
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Figure 7: Our algorithm detects 2 communities form Zachary karate club, which perfectly 
corresponds the actual division. 

B. College football network. The algorithm is also applied to College football net- 
work which was provided by Mark Newman. The network is a representation of the schedule 
of Division I games for the 2000 season: vertices in the network represent teams and edges 
represent regular-season games between the two teams they connect. What makes this 
network interesting is that it incorporates a known community structure. The teams are 
divided into 12 conferences. Games are more frequent between members of the same con- 
ference than between members of different conferences. The average shortest path length 
of the football network is 2.5, so we let the signaling time T = 3. When the F statistic 
achieve the maximum we get the best partition. We also use the accuracy index of detecting 
community algorithm [30J to measure the effect of our algorithm and find that it identifies 
the conference structure with a high degree of success. We detect 15 communities when F 
reach it's maximum (Figl8|) among which 2 communities just have 2 and 3 teams respectively 
and five communities were detected absolutely. The average accuracy is 0.78 which is little 
better than any of others (Tab l3.2p . 

4 Conclusion 

The investigation of community structures in complex networks is an important issue in 
many domains and disciplines. This problem is relevant to social tasks, biological inquiries, 
or technological problems. In this paper, we have introduced a method to detect communities 
based on the signaling process on networks. 

In a complex networks with n nodes, every node is viewed as a system which can be 
inspired. Each node sends its neighbors and itself signals and record the number of signals 
it receives at every time step. For each node of the network, we give it an initial unit quantity 
signal and other nodes have the signal of zero. Then after the inspiration of T steps on the 
network, the signal distribution of the nodes denoted by an n-dimensional vector can be 
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viewed as the influence of the source node to the whole network. The amount of signals be 
sent is equal to its present-owning signals. In complex networks, we can generally consider 
that the node always influence its community first then through the community influence 
the whole network. So naturally, compared with the nodes of other communities, nodes 
in the same community have similar influence toward the whole network. So through the 
signaling process, the network partition problem is transformed into the vectors clustering 
problem in Euclidian space. The clustering can be work out by Fuzzy C-means clustering 
method(FCM) with the help of F statistic. Moreover, our algorithm can also be generalized 
to weighted networks when we think the weighted connections can magnify or dwindle 
the signals linearly. Thus the method presented here can detect the optimal community 
structures in binary and weighted networks without any extra parameters. 

To solve the partition problem of complex networks, precision and accuracy of an algo- 
rithm are two standards for us to choose the method. So we make a comparison between 
our algorithm and other relatively mature ones such as EO, Potts and GN algorithms both 
in binary and weighted networks. Results for both ad hoc and real networks have proved 
that our algorithm is effective. One problem of our algorithm is that we haven't given a 
clear range of the most optimal steps T of inspiring. Actually, some other algorithms such 
as the random walks method [14^ exist also the similar problem. So we think it is an open 
problem and we will do some deep research on it in the future. 
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